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Abstract 

The problem of the possible creation of mixed hadron-quark-gluon matter, that 
can arise at nuclear or heavy-ion collisions, is addressed. It is shown that there can 
exist several different kinds of such a mixed matter. The main types of this matter 
can be classified onto macroscopic mixture, mesoscopic mixture, and microscopic 
mixture. Different types of these mixtures require principally different descrip- 
tions. Before comparing theoretical results with experiments, one has to analyze 
thermodynamic stability of all these mixed states, classifying them onto unstable, 
metastable, and stable. Only the most stable mixed state should be compared with 
experiment. Mixed states also need to be checked with regard to stratification insta- 
bility. In addition to the static stratification instability, there can happen dynamic 
instability occurring in a mixture of components moving with respect to each other. 
This effect, called counterflow instability, has also to be taken into account, since it 
can lead to the stratification of mixed matter. 
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1 Introduction 



At high temperatures and/or densities, hadronic matter is expected to undergo a transi- 
tion to quark-gluon plasma, where quarks and gluons are no longer confined inside hadrons 
but can propagate much further in extent than the typical sizes of hadrons. Such a de- 
confinement transition can happen under heavy-ion or nuclear collisions. It is assumed to 
exist in the early universe cosmology, since for a time on the order of the microsecond the 
temperature was high enough for the elementary degrees of freedom of QCD to be in a 
deconfined state. The quark-gluon plasma can also exist in the interior of compact stars. 

The peculiarities of the transition from hadronic matter to quark-gluon plasma, that 
is, of the deconfinement transition, have been the object of many discussions (see review 
articles [1-10]). From general arguments, it is impossible to infer the order of the QCD 
transition, whether it is 1-st order, 2-nd order, or crossover. Being based on a model 
consideration, the deconfinement was shown to be a gradual crossover [6,7]. At the present 
time, this result has been confirmed by numerical simulations of lattice QCD showing 
convincingly that deconfinement is really a crossover [11-15]. 

In order to be able to describe the states of matter and phase transitions in ther- 
modynamic terms, it is required that the matter be at least in quasi-equilibrium. The 
experimental lifetime of fireballs, formed under heavy-ion collisions, is of order t exp ~ 1CT 22 
s [16,17]. The local equilibration time of nuclear matter is t\ oc ~ 1CT 23 s [18,19]. Since 
Uoc "C t exp , equilibration is feasible and thermodynamic language is applicable to treating 
the fireball states. 

A plausible assumption is that in the process of the transformation of hadronic matter 
into quark-gluon plasma there can arise an intermediate state of matter representing a 
mixture of hadronic and quark-gluon states [20,21]. Note that the manifestation of quark 
degrees of freedom, resulting in the appearance of the Blokhintsev fluctons [22], Baldin 
cumulative effect [23], and in the formation of multi-quark clusters, has also been assumed 
to occur even at temperatures essentially lower than the deconfinement point [24-31]. 

However, the nature of the mixed hadron-quark-gluon state has not been well un- 
derstood. It is the aim of the present paper to explain that, actually, there can exist 
several kinds of such a mixed state, with the main three types that can be classified onto 
macroscopic, mesoscopic, and microscopic mixed states. These states have rather different 
properties and require essentially different theoretical description. 

In the paper, we use the system of units, where the Planck and Boltzmann constants 
are set to one. 

2 Macroscopic mixed state 

This type of mixed state would arise if the deconfinement transition would be of first 
order [32-34]. Then, at the phase-transition point, the system, say fireball of a linear 
size L, separates into macroscopic domains of size / corresponding to hadron phase and 
quark-gluon phase, so that 

a</~L, (1) 

where a is mean interparticle distance. The domains are macroscopic, being of order 
of the system size L. They also are called droplets or blobs, or bubbles [35-40]. Their 
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topology is similar to the droplets of nucleons arising in the low-density nuclear matter 
[41,42]. The domains of different phases correspond to different vacua [43-46], with the 
physical states of different domains being mutually orthogonal [47]. 

Strictly speaking, under first-order phase transition, mixed phase can occur only at 
the transition point, where two kinds of pure phases meet each other, pure hadron phase 
and pure quark-gluon phase. The qualitative behavior of the transition temperature T c 
as a function of barion density ps is shown in Fig. 1. Hadron matter consists of only 
hadrons, interacting with each other through hadron-hadron interactions [48,49]. Pure 
quark-gluon plasma is described by an equation of state for free quarks and gluons, with 
taking into account their interactions [50] by incorporating some non-perturbative effects 
[51-53]. 

The straightforward order parameters are the density of hadron matter, ph and the 
density of quark-gluon plasma, p q . In the hadronic matter 

Ph > , p q = , (2) 

while in the quark-gluon plasma 

Ph = 0, p q >0. (3) 

It is also possible to use as an order parameter the Wilson loop [54,55]. 

Each type of particles is characterized by barion number B>i and strangeness S^. For 
simplicity, the particles are assumed to be neutral. The chemical potential of the i-type 
particles is expressed as 

Pi = p B Bi + p s Si (4) 

through the barion, p B , and strangeness, ps chemical potentials. The barion and strangeness 
densities are given by the relations 

PB = J2B tPt = ^-, Ps = J2S tPt = ^-, (5) 

in which P is pressure and pi is the density of the i-type particles. The behavior of 
pressure, under first-order phase transition, is shown in Fig. 2. Since the grand potential 
is Q = —PV, where V is the system volume, the larger pressure corresponds to the lower 
grand potential. 

The transition temperature is defined by the equality of the pressures, 

P h (T c , p B ) = P q (T c , p B ) , (6) 

where, for simplicity, the strangeness density is fixed. This gives T c = T c (pb). Because 
from the left and the right of T c , the pressures are different, we have two barion densities, 
for hadrons and for plasma, 

pBh = = pBh(T, p B ) , p Bq = tt^- = pBq(T, p B ) , (7) 

Ops Ops 

which gives two barion potentials, psh and ps q that coincide at the transition tempera- 
ture: 

PBh{T c ,PBh) = PBq{T c ,PBq) ■ (8) 
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This defines T c = T c (p Bh ,p Bp ). 

The point of a first-order phase transition is the point of instability. Infinitesimally 
small fluctuations of temperature around T c will result in finite jumps between two dif- 
ferent barion densities in Eq. (J7|). So that the mixed phase at this point is unstable. 

One says that the mixed phase could exist not merely at the transition point, but also 
in a region around it. This is explained as being due to the Maxwell construction that is 
demonstrated in Fig. 3 for the pressure as a function of the reduced barion volume 

1 

v B = — • 

PB 

Here, the standard behavior of the pressure under a first-order phase transition [56] is 
corrected by replacing the part, corresponding to unstable and metastable states (shown 
by the dashed line), by the horizontal solid line between the barion volumes 

1 1 

VBh = , V Bq = • 

PBh PBq 

As a result of this construction, the phase diagram of Fig. 1 transforms into that 
of Fig. 4. The Maxwell construction for the pressure as a function of temperature is 
equivalent to the smoothing of the pressure, as is shown in Fig. 5. The mixed phase 
exists between the low, T n , and upper, T*, nucleation temperatures. 

However, as is evident from Fig. 3, on the coexistence line, one has 

dP 

= (v B h < V B < V B q) ■ (9) 

This implies that the isothermal compressibility 



1 / <9P x ~ 



vb \dv Bi 

is divergent everywhere in the region of the mixed phase existence: 

k t oo (T„ < T < T* n ) . (10) 

The divergence of the compressibility means instability, since infinitesimally weak pressure 
fluctuations would lead to the system explosion during the short explosion time 



texp r * J 
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That is, the fireball would explode even before it could equilibrate. 

Concluding, if deconfinement would be a first-order phase transition, then, formally, 
a mixed hadron-quark-gluon phase could arise around the transition point, however such 
a mixed state is strongly unstable and, in reality, cannot exist as an equilibrium phase. 
In addition, as QCD lattice simulations prove [11-15], deconfinement is not a first-order 
transition, but rather a crossover. 
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3 Mesoscopic mixed state 



There can exist another type of mixed state that can be called mesoscopic mixed hadron- 
quark-gluon matter. The term "mesoscopic" means that the typical size l g of the arising 
germs of one phase inside the other is between the mean interparticle distance and the 
system size: 

a</ 9 <L. (11) 

Below the deconfinement temperature, these are the germs of quark-gluon plasma sur- 
rounded by hadron matter. And above the transition temperature, these are the germs 
of hadron matter inside quark-gluon plasma. 

For the mesoscopic mixed state, pressure is uniquely defined and does not require the 
phase transition to be of the first order. Generally, it can be of any order, including the 
crossover type [19,57]. Mesoscopic mixed state can exist in a large temperature interval 
between the low and upper nucleation temperatures. 

The mesoscopic mixed state is basically different from the macroscopic one, exhibiting 
the following main features. 

(i) The germs of competing phases do not need to be in absolute equilibrium. They 
can have finite lifetime t g . But they are to be in quasi-equilibrium, such that the local 
equilibration time be essentially shorter than their lifetime: 

tl0C < tg . (12) 

(ii) The spatial distribution of germs at a snapshot is random. They form no ordered 
spatial structure, such as domains. 

(iii) The spatial distribution of germs is also random with respect to repeated experi- 
ments. 

(iv) The typical size of the germs is mesoscopic in the sense of Eq. ffTTl) . 

(v) The germ geometry is of multiscale nature. Their shapes are not regular, but are 
rather ramified. And the sizes l g lie in a dense interval [Z™ n , l g nax ], such that 

a < Z™ n < l™ ax < L . 

The description of the mesoscopic mixed state has to take into account the generic 
random nature of the spatial germ distribution. The main ideas of the theory are as 
follows [19,57,58]. We keep in mind a mixture of two phases, e.g., hadron matter and 
quark-gluon plasma. 

At a snapshot, the system volume is divided onto the volumes of different phases, 

V = Vl |jv 2 , 

separated by the Gibbs equimolecular separating surface, for which extensive observable 
quantities are additive. This concerns as well the number of particles in each phase and 
the related volumes: 

N = Ni + N 2 , V = Vy + V 2 , (13) 

with V v = mesV^. Mathematically, the separation is characterized by the manifold indi- 
cator functions 
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where r is a spatial variable and p — 1,2 enumerates the phases. 

At a snapshot, the mixture needs to be described by a representative statistical en- 
semble {H,p(Q}, where W is the space of microstates and />(£) is a statistical operator 
[58]. The space of microstates is given by the fiber space 



1-2 , 



(15) 



with the fiber bases T-L v being weighted Hilbert spaces. The statistical operator is nor- 
maized as 

(16) 



Tr J Vi = 1 



by taking the trace over the quantum degrees of freedom and averaging over the random 
germ spatial configurations defined through the functional integral over the manifold 
indicator functions (I14p . 

To construct a representative ensemble, one defines the internal energy 



and all constraining quantities 



e = Tr J mm) V£ 
a = Tr J mciio ■ 



(17) 



;is) 



required for the unique description of the system. The statistical operator is found from 
the principle of minimal information, by minimizing the information functional 



/[/)(£)] = Tr J p(0 lnp(0 Vi + 



+ A 



Tr / p(0 Vi-1 



Tr / p(£)H(Q Vi-E 



+ 



Tr / p(0Ci(0 Vi - C t 



(19) 



in which Ao, 0, and Aj are Lagrange multipliers. 
The minimization yields the statistical operator 



m = % exp{-/3i?(0} 



with the grand Hamiltonian 



i 

where /i, = —\T. The partition function is 



(20) 



(21) 



Z = Tr / exp{-/3#(0} V£ 
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and (3 = 1/T is inverse temperature. 

Let us introduce the effective Hamiltonian H defined by the equality 



exp{-/3# (0} Vi = exp(-/3F) . (22) 
After this, the partition function reduces to the form 



Z = Tre 



-fSH 



containing only the trace over quantum degrees of freedom. 

The geometric weights of each phase are given by the expressions 



w„ = J £,(r) Vi, (23) 

satisfying the normalization condition 

w x + w 2 = 1 . (24) 
This provides the minimum for the grand potential 

n = -TlnTre _/ji? , 

that can be found from the conditions 

dQ d 2 Q . . 

— - = — > , 25 
aw v owf, 

taking into account normalization (l24p . The phase weights f[2^|) play the role of additional 
order parameters characterizing the mixed state [19,57,60,61]. 

The mesoscopic mixed state is stable, with deconfmement being rather a sharp crossover. 



4 Microscopic mixed state 

The third type of mixture is termed microscopic because hadrons are uniformly inter- 
mixed with quark-gluon plasma, without forming either germs or droplets. Such a mixed 
state can be treated by the theory of clustering matter [6,7], considering hadrons as quark 
clusters. Each kind of clusters, enumerated by the index i, is characterized by the bar- 
ion number Bi, strangeness Si, and compositeness Zi. The latter shows the number of 
quarks forming a cluster of that type. For instance, the quark compositeness is 1, meson 
compositeness is 2, and the nucleon compositeness is 3. 

The space of microstates for the mixture is the tensor product 

M = ®Fi, (26) 

i 

in which Ti is the Fock space for the i-clusters. 
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The density of i-clusters is 



Pi = d J ni(k) , (27) 



where Q is a degeneracy factor and nj(k) is a momentum distribution. The total mean 
quark density is 

p = J>Pi. ( 28 ) 

j 

The cluster weights are defined by the ratio 

Wi=^, (29) 
P 

which gives Wi = Wi(ps, Ps,T). By definition, one has 

< Wi < 1 , W{ — 1 . 

i 

The Hamiltonian of a microscopic mixture, generally, has the form 

^ = E^ + ^E^-' ( 3 °) 

i " i^j 

in which the first term is the sum of the channel Hamiltonians and the second term 
corresponds to cluster interactions. Modeling the Hamiltonian, one often assumes its 
dependence on density and/or temperature. For example, the effective particle spectra 
are often defined as functions of temperature [62]. Therefore, in the definition of the 
grand Hamiltonian, 

H=H-^fHNi + CV , (31) 

i 

one has to include the term CV guaranteeing statistical correctness for the approach. To 
this end, it is necessary to require the validity of the conditions 



\w °' wr ' (32) 



The latter reduce to the equations 

QH\ Fir, 1 / r)H\ 

(33) 



dC 1 /9H\ DC 1 IdH 

dp~~~V\'dp~/ dT~~V\dT 

defining C = C({ Pl },T). 

Only under the conditions of the statistical correctness (!3~2]) . the theory becomes self- 
consistent and satisfies all thermodynamic relations: 

P = -V = ~W E= V {H)=T dT~ P + WB + t*sPs, 

dP 1 1 - dP 

S= df = f (yE + P ~ PbPb ~ ^ sPs > ' Pi= V ( Ni ' = djT- 



p B = > BiPi = - — , p s = > biPi = - — 
^ op B ^ dps 



It is a common mistake, widely spread in literature, when the authors forget about sta- 
tistical correctness, because of which the obtained results cannot be reliable. 

Taking into account cluster interactions may seem to be a problem, since there can 
exist various quark clusters, whose interactions are not known. This obstacle can be 
avoided in the following way. Let us consider the reaction of fusion of two clusters, say 
a cluster a and cluster b, into one cluster i, with all compositeness numbers larger than 
one, so that there is the conservation of compositeness, 

Za ~\~ %b i 

and the conservation of mass, 

m a + m b + $ a{ , = rrii , 

where <3> a {, is the interaction energy of two clusters. For the same fusion, in the presence 
of a third cluster j, the mass conservation reads as 

m a + m b + rrij + $ a6 + $ aj + $ 6i = m, + rrij + ^ . 

From these relations, it follows the potential scaling law 

(34) 



Z{Zj Z a Z\y 

This law allows us to express all cluster interactions through one known interaction, e.g., 
through the nucleon-nucleon interaction, 

= *ZL<f> NN , (35) 

which is well known [49]. 

The microscopic hadron-quark-gluon mixture is stable, with deconfinement being a 
sharp crossover [6,7], in good agreement with the QCD lattice simulations [11-15]. In the 
case of a crossover, the deconfinement temperature can be defined as the point where the 
derivatives of observables have a maximum, which gives about 170 MeV. Of course, con- 
sidering different observables can result in slightly different deconfinement temperatures, 
which is the common situation for crossovers, where the crossover temperature is defined 
conditionally. Numerical simulations [63,64] show that pion clusters survive till around 
2T C . 



5 Static and dynamic stability 

One more problem that arises in considering the coexistence of clusters of different types 
is the possibility of their spatial stratification, when the clusters, previously uniformly 
mixed, separate in space into domains containing only one kind of clusters. Below, we 
illustrate this problem by considering a two-component mixture of clusters. 
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Let the total number of clusters be N = Ni + N2, existing in the volume V = V\ + Vi- 
The system can form two kinds of mixture. One situation corresponds to a microscopic 
mixture, with all clusters being uniformly intermixed in the space. And the other case is 
when the clusters of each type are spatially separated into different domains, thus forming 
a macroscopic mixed state. The microscopic mixture is more thermodynamically stable 
when its free energy F mix is lower than the free energy F sep of the separated state of the 
macroscopic mixture, 

F m ix ^ F se p . (36) 

Calculating the free energy in the correlated mean-field approximation, we use the nota- 
tion for the mean interaction intensity 

% = J ^(r)^(r) dr , (37) 

in which V^-(r) is a vacuum cluster interaction and <7y(r) is the pair correlation function. 
Then from Eq. (IHtjj) . we find the condition for the stability of the microscopic mixture 

/ TV 

$12 < V*n*2 2 + AS mix , (38) 

iVliV2 

where AS m i X is the entropy of mixing, which can be written as 

iVi N 9 
AS mix = -N! In - N 2 In . (39) 

In the case of validity of the potential scaling ( l34|) . the stability condition fl38|) reduces 
to the trivial requirement that the entropy of mixing (139]) be positive, which is certainly 
true. Hence, under the validity of the potential scaling, the microscopic mixture is always 
more stable and there is no stratification. 

The stability condition ( 138)) is derived for an equilibrium situation by comparing the 
thermodynamic potentials of the microscopic mixture and the separated stratified state. 
In that sense, it is a static stability condition. But there is also a dynamic stability con- 
dition requiring that the spectrum of elementary excitations be positive [65]. Analyzing 
the dynamic stability, we take into account that the components can move with respect 
to each other with the velocities vi and V2. Such a relative motion can be due to the fact 
that the fireball has been formed as a result of two colliding heavy ions or nuclei. 

Studying the spectrum of collective excitations of a microscopic mixture in the random- 
phase approximation, we find that the spectrum is positive, provided that the relative 
velocity v = v 2 — Vi does not exceed by the magnitude v = |v| the critical value 



''- ($ n $ 22 - $? 2 ) . (40) 



m 2 $ 
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If v < v c , the microscopic mixture is stable. But if v > v c , the mixture stratifies into 
macroscopic domains containing different sorts of clusters each. The dynamic instability, 
leading to the stratification, caused by the mutual motion of components, is called the 
counterflow instability. 

In conclusion, we have explained that there are three types of mixed systems, macro- 
scopic, mesoscopic, and microscopic. Each kind of these mixed states is very different 
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from others, enjoying quite different physical properties and needing principally different 
theoretical description. 

If deconfinement would be a first-order phase transition, there could arise the macro- 
scopic mixed state, where hadron and quark-gluon phases would be located in separate 
macroscopic spatial domains. However, such a state is not stable and would disappear 
even before a fireball would equilibrate. In addition, lattice QCD simulations demonstrate 
that deconfinement is not a first-order transition, but a crossover. Hence, the macroscopic 
mixed state has no chance to exist. Therefore the naive picture, when one compares the 
mixed hadron-quark-gluon phase with a boiling water containing gas bubbles, has nothing 
to do with QCD. Theoretical predictions, based on the macroscopic mixed model, cannot 
be confronted with experiment. 

The real quark-hadron mixed state can be either mesoscopic or microscopic. These 
states can be stable, with deconfinement being rather a sharp crossover. 

Studying a multicomponent mixture, it is necessary to check it with respect to the 
stratification instability. The components, moving through each other, can also exhibit the 
counterflow instability. All these effects need to be carefully analyzed before comparing 
theoretical predictions with experimental observations. 
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Pb 

Figure 1: First-order phase-transition temperature T c as a function of barion density. 




Figure 2: Pressures of the hadron matter and quark-gluon phase as functions of temper- 
ature. 
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Figure 3: Maxwell construction for the pressure as a function of the reduced barion 
volume. 



T 




Figure 4: Mixed hadron-quark-gluon phase formally appearing around a first-order phase 
transition result of the Maxwell construction. 




Figure 5: Smoothed pressure (dashed line) as a function of temperature, corresponding 
to the Maxwell construction. 
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